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A molecular Rayleigh scattering technique is developed to measure dynamic gas 
temperature, velocity, and density in unseeded turbulent flows at sampling rates up to 10 
kHz. A high power CW laser beam is focused at a point in a heated air jet plume and 
Rayleigh scattered light is collected and spectrally resolved. The spectrum of the light, 
which contains information about the temperature, velocity, and density of the flow, is 
analyzed using a Fabry-Perot interferometer. The circular interference fringe pattern is 
divided into four concentric regions and sampled at 1 and 10 kHz using photon counting 
electronics. Monitoring the relative change in intensity within each region allows for 
measurement of gas temperature and velocity. Independently monitoring the total scattered 
light intensity provides a measure of gas density. Power spectral density calculations of 
temperature, velocity, and density fluctuations, as well as mean and fluctuating quantities 
are demonstrated for various radial locations in the jet flow at a fixed axial distance from 
the jet exit plane. Results are compared with constant current anemometry and pitot probe 
measurements at the same locations. 
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Nomenclature 

most probable molecular speed (ms -1 ) 
speed of light (= 2.998 x 10 8 m s" 1 ) 

Fabry-Perot mirror spacing (m) 
incident electric field vector (N C" 1 ) 

Fabry-Perot instrument function contrast 
frequency of scattered light (s' 1 ) 
frequency of incident laser light (s' 1 ) 
collimating lens focal length (m) 
fringe forming lens focal length (m) 
sampling rate (s' 1 ) 

Planck’s constant (= 6.626 x 10" 34 N m s) 
Fabry-Perot instrument function 
interaction wave vector (m -1 ) 
magnitude of K (m -1 ) 
incident light wave vector (m -1 ) 
scattered light wave vector (m -1 ) 
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length of sub-records in spectral calculations 
probe volume length (m) 

mean free path between molecular collisions (m) 
molecular mass (kg) 
effective finesse 

photoelectron counts at q th annular region 

Rayleigh scattered photoelectrons collected into solid angle Q 

molecular number density (nf 3 ) 

power of incident laser beam (N m s" 1 ) 

static pressure (N m" 2 ) 

annular region number 

radial position in image plane (m) 

radial position in the jet measurement plane (m) 

fiber face image radius (m) 

inner radius of q th annular region (m) 

fringe radius associated with reference laser light (m) 

Rayleigh scattering spectrum 
static temperature (K) 
velocity vector (ms -1 ) 
measured velocity component (ms -1 ) 

frequency shift non-dimensionalized by the frequency shift related to the local speed of sound 
y-parameter 

angle between E and k s vectors (rad) 

Fisher information matrix 

width of annular region (m) 

integration time (s) 

optical system efficiency 

dynamic viscosity (N s m" 2 ) 

angle between light ray and optical axis (rad) 

angle between reference light ray and optical axis (rad) 

Boltzmann’s constant (= 1.381 x 10" 23 J K" 1 ) 
wavelength corresponding to wave vector K (m) 
illumination wavelength (m) 
refractive index 
gas density (kg m" 3 ) 

differential scattering cross-section (m 2 sr" 1 ) 
scattering angle (rad) 

phase change between successive reflections 
solid angle (sr) 


I. Introduction 

T here is currently no high frequency response, non-intrusive temperature measurement technique for use in 
turbulent flow studies. There are also no techniques that can provide non-intrusive dynamic measurement of 
gas temperature, velocity, and density simultaneously. Conventional intrusive measurement devices such as 
resistance wires, pressure probes, and thermocouples are usually limited in spatial and temporal response, disturb the 
flow under study, and can be damaged by high pressure or temperature. Therefore, a non-intrusive molecular 
Rayleigh scattering technique is developed to measure dynamic gas temperature, velocity, and density in unseeded 
turbulent flows at sampling rates up to 10 kHz. 

Various optical techniques are available for temperature, molecular number density, and velocity measurements. 
Eckbreth 1 provides the details of several techniques, such as Coherent Anti-Stokes Raman Spectroscopy (CARS), 
Spontaneous Raman and Rayleigh scattering, and Laser-Induced Fluorescence (LIF), with an emphasis on their use 
in combustion applications. Other techniques that have been used for flow measurements in various environments 
include Laser-Induced Thermal Acoustics (LITA) 2 ’ 3 , also known as Transient Grating Spectroscopy (TGS) 4 , Laser- 
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Doppler Velocimetry (LDV) 5 " 8 , filtered Rayleigh scattering 9 , Particle Image Velocimetry (PIV) 10 , and focusing 
schlieren deflectometry 11 . However, none of these techniques can provide high frequency response measurement of 
all three flow properties simultaneously. 

The temporal response of most of the optical measurement techniques mentioned here is typically limited by the 
repetition rate of pulsed lasers, usually on the order of 10 Hz; however one can potentially take many snap shots (1-5 
ns temporal resolution) to obtain turbulence measurements. Signal strength of the scattering process can also limit 
the temporal response. Rayleigh scattering has several advantages over other light scattering techniques. The 
Rayleigh scattering cross-section is three orders of magnitude greater than the Raman scattering cross-section, and 
Rayleigh scattering is generally much simpler to implement than techniques that require multiple lasers and have 
other complicating factors. 

Previous works using molecular Rayleigh scattering to make temperature, velocity, and number density 
measurements in harsh environments have been reported 12 " 15 . The current work is an extension of a previously 
developed technique using Rayleigh scattering to obtain dynamic density and velocity measurements in supersonic 
free jets 16,17 and dynamic density and temperature measurements in heated jets 18 . The Rayleigh scattering technique 
presented is unique in that it provides simultaneous temperature, velocity, and density fluctuation measurements. A 
high power CW laser beam is focused at a point in a jet plume and Rayleigh scattered light is collected and 
spectrally resolved. The spectrum of the light contains information about the temperature, velocity, and density of 
the flow. A planar mirror Fabry-Perot interferometer (FPI) is used to analyze the spectrum of the scattered light, and 
photomultiplier tubes (PMTs) record the fringe intensity pattern at high sampling rates enabling dynamic 
measurement of the gas flow properties. 

The ability to obtain high frequency response measurement of multiple properties simultaneously is a valuable 
tool for compressible, turbulent flow research. These data are particularly useful to aeroacoustics researchers who 
are interested in correlating flow property fluctuations with far field acoustic fluctuations. Knowledge of velocity- 
temperature fluctuation correlations is also needed to improve computational fluid dynamic (CFD) models of 
compressible turbulent flows. 


II. Theory 


A. Rayleigh scattering and spectral analysis 

Molecular Rayleigh scattering is the result of elastic light scattering from gas molecules. When light from a 
single frequency laser beam passes through a gas, the scattered light is shifted in frequency by the Doppler effect 
due to the thermal as well as the bulk motion of the molecules. The frequency spectrum of the scattered light 
contains information about the gas density, bulk velocity, and temperature. Figure 1 shows a Rayleigh scattering 
spectrum containing the narrow laser line and a typical Rayleigh spectral peak to illustrate how the flow property 
measurements are obtained from the spectral information. If the gas composition is fixed, the total intensity of the 
Rayleigh scattered light is directly proportional to the gas density. The frequency shift between the laser peak and 
the Rayleigh peak is proportional to the bulk flow velocity. The width of the spectrum is related to the gas 
temperature. 




Figure 1. Rayleigh scattering spectrum. 


Figure 2. Light scattering from a moving particle. 
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The shape of the spectrum is dependent on gas pressure, temperature and the scattering angle 19 . A non- 


dimensional parameter y, which represents the ratio of the wavelength of the scattering grating 
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mean free path between molecular collisions [l m = 2^-j, is used to establish spectral shape regimes: 
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The interaction wave vector, which defines the direction of the velocity component being measured, is the 
bisector of the incident and scattered light wave vectors (Fig. 2). The interaction wave vector and its magnitude are 
given by: 


K - k s -k 0 
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The geometry of the optical arrangement in an experiment can be designed such that the desired component of the 
velocity is measured: 


u k = 


K u 
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( 5 ) 


Three spectral shape regimes are defined for typical 90° scattering 19 " 21 . For low density gases where y « 1, the 
Rayleigh spectrum is accurately modeled by a Gaussian function and the gas is said to be in the Knudsen or 
collisionless regime. For higher density gases where y » 1, the Rayleigh spectrum broadens and eventually 
develops side-lobes known as Brillouin peaks. This is known as the hydrodynamic or continuum regime where 
molecular collisions are the dominant process. Finally, for 0.2 <y < 2, as is the case in the present experiments, the 
gas is in a transition region between the collisionless and hydrodynamic regimes, and a kinetic theory model is 
required. The shape regimes are illustrated in Fig. 3, where the normalized frequency is defined as: 


2 n(f-h) 

f Vn 


( 6 ) 


A kinetic theory model developed by G. Tenti 19,22 provides a description of Rayleigh-Brillouin scattering from 
molecular gases in all density regimes. The TENTI S6 spectrum model 19,22 was used to generate the information 
displayed in Figs. 1 and 3, and is incorporated in the model function used in least squares analysis of the 
experimental data. 

The spectrum of the Rayleigh scattered light is analyzed using a planar mirror FPI (Fig. 4) operated in the static 
imaging mode 23 . The fringe intensity pattern is a function of the Rayleigh spectrum and the Fabry-Perot instrument 
function. The Fabry-Perot instrument function is: 
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Figure 3. Rayleigh scattering spectrum for Figure 4. Fabry-Perot Interferometer, 

various y-parameters. 


In our experiment, Rayleigh scattered light from a defined probe volume is collected into a multimode optical 
fiber. The fiber directs the light to the FPI and a lens at the interferometer output focuses the interference fringe 
pattern at the image plane. The total expected number of Rayleigh photoelectron counts without the FPI in the 
optical path can be expressed as: 


( n r) 


£ P 0 n L x X Cl At f dcr ^ 
he \ dCl ) 


sin 2 /? 


( 11 ) 


where the overall system efficiency £ includes detector quantum efficiency and other losses. When the 
interferometer is placed in the light path, the model function for the amount of energy collected from the q th annular 
region of the interference pattern, which has been dissected into one circular and three annular regions, in terms of 
photoelectron counts, can be expressed as follows: 



(Nj?) rq+Ar 00 

X 2 7 | \S R (x f )I FP (x f ,r)dx f dr 

ft fnax r q -oo 


( 12 ) 


where the Rayleigh spectrum is evaluated using the TENTI S6 model. 

B. Lower bound for measurement uncertainty 

The lower bound on the uncertainty in temperature, velocity, and density measurements using Rayleigh 
scattering is set by the photon statistical (shot) noise. Estimates of the measurement uncertainty in the unknown 
parameters for this technique were obtained by Cramer-Rao lower bound analysis 24 . For a measurement that is a 
function of a set of unknown parameters, a i , the variance of the estimates of the parameters is: 
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(13) 


where the Fisher information matrix for Poisson statistics is given by: 



(14) 


Using the model function developed for^/V^ , the uncertainty or standard deviation cr(cr z ) in each unknown 
parameter may be determined by calculation and inversion of the Fisher information matrix. 

If a Gaussian spectrum is assumed, the lower bounds for uncertainties in temperature T 6 7 , velocity u k , and density 
p for an ideal instrument are: 
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For air flow at T = 360 K, u k = 9 m/s, and standard pressure, the lower bound uncertainty estimates for a 
sampling rate of 10 kHz, assuming an ideal instrument, are: 

-^2=3.6% cr(u k ) = 8 m/s ^2=5.4% 

T P 

For any real instrument, these measurement uncertainties will be higher. The actual Cramer-Rao lower bound 
measurement uncertainties were calculated, given the Fabry-Perot interferometer, PMTs, and overall efficiency of 
the current system. The TENTI S6 spectral model 19,22 was used to calculate the Rayleigh scattering spectrum in this 
analysis. Figure 5 shows the uncertainty in T and u k over a range of fringe radii (r R ) related to the incident 
(reference) laser light. This analysis was used to determine the reference fringe radius resulting in minimum 
uncertainty levels in both T and u k . A r R value of 6.5 mm was chosen for the experiments. Figure 6 gives the 
uncertainty in T and p for temperatures ranging from 300 to 700 K, with u k = 9 m/s, r R = 6.5 mm, and sampling rate 
of 10 kHz. The measurement uncertainty of all parameters increases as temperature increases. This is attributable to 
the lower gas density at higher temperatures, resulting in less scattering molecules, and hence fewer scattered 
photons. Although the uncertainty in the parameters is rather high for instantaneous measurements, long data records 
allow for calculation of higher accuracy statistical quantities such as power spectra and mean square fluctuations. As 
sampling rate decreases, so will uncertainties, since the uncertainty is inversely related to the total number of photon 
counts, as shown in Eqs. 15. Experimental data was acquired at sampling rates of 1 and 10 kHz. The lower sampling 

rate increased photon counts by a factor of 10, thereby improving uncertainty levels by a factor of VlO over the 
higher sampling rate. 
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Figure 5. Uncertainty in velocity and temperature as 
a function of reference fringe radius; u k = 9 m/s, 
J=360K. 


Figure 6. Uncertainty in density and 
temperature as a function of temperature; u k = 9 
m/s, = 6.5 mm. 


III. Experiment 

Flow measurements were performed downstream of a heated air stream issuing from an 8 mm diameter straight 
tube surrounded by a 50 mm diameter 2 m/s room temperature co-flow at the NASA Glenn Research Center using 
the described Rayleigh scattering measurement technique. The core flow was heated via electrically heated coils 
inside the tube, which generated turbulent mixing in the exiting air stream. The heated air jet was operated over a 
static temperature range of 295 to 700 K with a mean centerline velocity ranging from 8 to 12 m/s measured at an 
axial distance of 32 mm (4 jet diameters) from the tube exit. The Reynolds number ranged from 1500 - 3000. A PID 
feedback control unit was used to adjust and maintain the jet total temperature to within +/- 1.5 K of the set-point 
temperature. Since the measurement technique relies on having particulate free gas flows, a series of filters were 
placed in line with the air plumbing to remove dust, oil and water from both air supplies. Figure 7a shows the layout 
of the optics around the jet, which were used to collect Rayleigh scattered light from gas molecules in the flow. The 
jet was mounted such that the main flow direction was parallel to the table. A 5W, 532 nm wavelength, single- 
frequency, Nd: Vanadate CW laser provided incident light for the system. The laser beam was focused with a 250 
mm focal length lens (LI) to a 70 pm 1/e 2 diameter at the probe volume positioned 32 mm downstream of the jet 
exit. The beam was oriented in the horizontal plane, at a 45° angle to the primary flow direction. Light was collected 
at a 90° scattering angle, collimated by a f/4 200 mm focal length lens (L2), and focused by a 100 mm focal length 
lens (L3) onto a 0.55 mm diameter multimode optical fiber. Since the pair of lenses provided 1:2 imaging, the length 
of the probe volume was 1.1 mm. The incident and scattering wave vectors were arranged such that the axial 
component of the jet velocity was measured, as shown in Fig. 7a. The jet was mounted on vertical and horizontal 
translation stages so that the 1.1 mm probe volume could be positioned anywhere in the jet plume. 

The fiber was routed to a separate optical table where spectral analysis and detection took place, shown in Fig. 
7b. The light exiting the fiber was collimated by an 80 mm focal length lens (L5) and was directed through the 
planar mirror FPL The FPI had 70 mm diameter mirrors with 80% reflectivity, 8.7 GHz free spectral range (FSR), 
and reflective finesse of approximately 16. The Fabry-Perot is an extremely sensitive instrument; even the smallest 
vibrations or temperature changes can cause the mirrors to drift out of parallel alignment, resulting in increased 
uncertainty. Therefore, a stabilization system was utilized to maintain parallelism of the mirrors during testing. 
Between Rayleigh scattering measurements, a mirror and diffuser (Fig. 7a) were placed in the beam path by linear 
actuators to direct some of the incident laser beam into the optical fiber and through the FPL A set of reflecting 
prisms mounted on a linear actuator (Fig. 7b) was positioned in the optical path at the output of the interferometer to 
direct the light from three regions of the interferometer mirrors to a video camera. Live video of the three fringe 
images was used in a feedback control loop to adjust the mirror positions using piezoelectric transducers until the 
three fringes were equal in diameter. This system was also used to set the fringe diameter of the incident reference 
light to 13 mm, as determined from the uncertainty analysis. Further details regarding the mirror stabilization system 
for the FPI can be found in Ref. 17. 

When flow measurements were collected the prisms, mirror, and diffuser were removed from the optical path 
and the light exiting the FPI was focused by a fringe forming lens (L6) having an effective focal length of 2584 mm, 
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Figure 7. Diagram of experiment setup: (a) laser and collection optics; (b) spectral analysis and detection 
equipment. 

which provided an 18 mm diameter image of the fiber face. The light was sampled via a series of various diameter 
45° elliptical mirrors (Fig. 7b). Mirror M8 sent the central 5 mm diameter section of the fringe image to PMT 1. 
Lens L8 refocused the remaining portion of the image with a magnification factor of 0.87 at mirror M9 which sent 
the central 1 1.5 mm diameter section of the original image to PMT 2. A pair of lenses (L10 and LI 1) refocused and 
de-magnified the image again so that mirror M10 removed a circular section equivalent to the central 16.1 mm 
diameter region of the original image, and sent this light to PMT 3. The remaining outermost annular region of the 
image was focused at PMT 4. Figures 8 and 9 illustrate the individual areas of the fringe image which are directed 
towards each PMT. The PMTs were operated in the photon counting mode to acquire fringe intensity data. The 
photoelectron pulses from the five PMTs were amplified (Gain = 5) and sent to photon counters, which output NIM 
level pulses. Constant fraction discriminators converted the NIM pulses to 10 ns wide TTL level pulses that were 
counted via counter-timer boards. Typical photoelectron count rates for this work were on the order of 1 MHz. A 
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Figure 8. Dissection of Fabry-Perot fringe Figure 9. Radial slice of a fringe pattern showing the 

pattern into four annular regions. annular regions detected by the four PMTs. The amount of 

light in each region changes as the velocity and 
temperature fluctuate. 

PC-based data acquisition system was used to record the signals from the PMTs. Each signal channel was digitized 
at either 1 or 10 kHz sampling rate for 30 seconds, resulting in 30,000 or 300,000 samples per channel. The intensity 
information from the Fabry-Perot fringe sampling was used to obtain temperature and velocity measurements. The 
density measurement was much simpler, requiring only an overall intensity measurement. A beamsplitter located 
just in front of the FPI input (Fig. 7b) directed approximately 10% of the incoming Rayleigh scattered light to a lens 
(LI 4), which focused the light at PMT 5 for density measurement. 

The temperature, density, and velocity measurements at the probe location were verified using additional 
instrumentation. A fine-wire (0.25 mm diameter) open-bead type K thermocouple provided jet static temperature 
data. A pressure gauge measured the ambient (static) pressure. Gas density information was derived using the ideal 
gas law. The jet velocity profile was evaluated from total pressure measurements with a pitot probe. A constant 
current anemometry (CCA) system provided static temperature fluctuation measurements. The system consisted of a 
1.3 pm diameter, 0.9 mm long platinum resistance wire operated in a constant current configuration. The frequency 
response of the CCA system, based on the 3 dB point, was measured to be approximately 1 kHz. 


IV. Results 


A. Calibration and technique verification 

The technique was calibrated over a jet static temperature range of 295 - 700 K and velocity range of 8 - 12 m/s. 
Thermocouple and pitot probe measurements were used to determine the mean temperature and velocity of the gas 
at the center of the jet plume 32 mm downstream of the tube exit. Since thermocouples measure the temperature of 
the metal probe bead rather than the gas temperature, there may be a difference between the actual gas temperature 
and the thermocouple measurement. These differences result from conduction losses down the thermocouple wire 
and radiation losses to the surroundings. The error in thermocouple-measured gas temperature was estimated to be 
less than ±3 K for the temperatures studied here. Gas density was derived from thermocouple and pressure probe 
measurements. The photoelectron counts recorded by PMT 5 are linearly related to density as shown by the 
calibration curve in Fig. 10 for data acquired at two sampling rates (1 and 10 kHz). The photon counts were adjusted 
by a small amount to account for stray laser light and provide an intercept that indicates nearly zero photon counts at 
zero density. 
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Figure 10. Density as a function of PMT 5 photon 
counts. 


The axial velocity component and temperature were measured simultaneously by analyzing the light from four 
annular regions of the interference fringe pattern. A model of the fringe pattern, including the Rayleigh scattering 
spectrum and FPI instrument function, was used in a least squares fitting routine to estimate unknown parameters 
from the intensity data. Mean temperature, velocity, and pressure measured by physical probes were used as known 
values in the least squares analysis of the mean photon count data at the calibration points, while several unknown 
factors were evaluated over the range of temperatures and velocities. These unknown factors included system 
detection efficiencies, effective finesse of the FPI, and fringe forming lens focal length. Once these factors were 
calibrated, the mean PMT counts acquired at sampling rates of 1 and 10 kHz were analyzed to evaluate the mean 
temperature and velocity at each data point, with the constraint that the pressure was equal to the ambient value 
measured by the pressure probe. Figures 11 and 12 respectively show thermocouple and pitot probe measurements 
compared with temperature and velocity measurements from the Rayleigh technique over the range of calibration 
temperatures. The temperature measurements agree extremely well for the 1 kHz sampling rate, but start to deviate 
at the high end of the temperature range for the 10 kHz sampling rate due to extremely low photon counts. The 
velocity measurements appear quite poor at a first glance at Fig. 12. Keep in mind that the measurement uncertainty 
is on the order 10 m/s, which is approximately the same magnitude as the velocity being measured, so the relative 
uncertainty is very high. Error bars have been shown to represent the standard deviation of the instantaneous 
measurements. At higher velocities (> 100 m/s), 10 m/s uncertainty is not as significant. Although the velocity 
uncertainty is relatively high for the low velocities studied here, the results are shown to verify that simultaneous 
measurement of temperature and velocity is feasible. 
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Figure 11. Mean temperature measurements from 
Rayleigh technique compared with thermocouple 
measurements. 


Figure 12. Mean velocity measurements from Rayleigh 
technique compared with pitot probe measurements. 
The error bars represent the standard deviation of the 
instantaneous Rayleigh-measured velocity values. 
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Using the linear equation relating PMT 5 photon counts to density (Fig. 10), and least squares analysis of the 
photon counts from PMTs 1-4, time histories of instantaneous density, temperature, and velocity were evaluated at 
each data point. The power spectra and mean square fluctuations were calculated from this information. Figures 13 
and 14 respectively show the first 5 seconds of the time history temperature results from the CCA measurements and 
the Rayleigh measurements acquired at 1 kHz sampling rate. The 2 second cycles in the data were due to the 
temperature controller time setting of one second on, one second off. The Rayleigh data demonstrates higher noise 
levels due to shot noise in the PMT signals; however the results verify that the Rayleigh temperature measurements 
sufficiently reproduce the CCA results. 
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Figure 13. Time history of instantaneous 
CCA temperature measurements. 
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Figure 14. Time history of instantaneous 
Rayleigh temperature measurements. 


B. Power Spectral Density (PSD) calculations 

The Rayleigh probe volume and the CCA probe were scanned vertically across the jet centerline at an axial 
station of 32 mm. Scans were performed at centerline jet conditions of 360 K and 9 m/s. These radial scans provided 
measurements in the jet core as well as in the turbulent mixing layer formed between hot and cold air streams. 
Because of the shot noise in the PMT signals, it was necessary to use relatively long data records and calculate 
power spectral densities using a technique such as the Welch method of modified periodograms 25 . Due to extensive 
data processing time (approximately 700 samples/min), only the first 6 seconds of the 10 kHz data were analyzed. In 
the Welch method, a data record sampled at a rate f s = 10 kHz (1 kHz) for a total time of 6 sec (30 sec) was 
subdivided into smaller records of length L = 1024 (128) samples, which were overlapped by 50%. The modified 
periodograms of each sub-record were calculated using a data window. These individual periodograms were then 
averaged to obtain the estimate of the power spectrum. The frequency resolution of the resulting spectrum is f/L = 
9.77 Hz (7.81 Hz). Overlapping the segments by 50% provided a near maximum reduction in the variance in the 
spectral estimate. The resulting power spectral densities provide fluctuation information out to 5 kHz (0.5 kHz). The 
sum of all points in the PSD is equivalent to the mean square fluctuations. 

Figures 15-17 show temperature PSD plots for three radial measurement locations: centerline (r 7 = 0), shear layer 
(rj = 6 mm), and a region outside the mixing layer (r 7 = 9 mm), calculated from the 1 kHz CCA data and the 1 and 
10 kHz Rayleigh data. The shot noise results in a constant noise floor in the spectra. This noise floor was apparent in 
the 10 kHz data at frequencies above 2.5 kHz where the spectrum is very flat (no fluctuations). The noise floor could 
not be estimated from the 1 kHz data since the power spectrum does not provide information at higher frequencies 
where it is safe to say that there is no contribution from true temperature fluctuations. Therefore, the shot noise 
contribution has only been subtracted from the 10 kHz spectra. Comparison between the temperature PSD results 
obtained from the Rayleigh data and the CCA data in the 0-0.5 kHz range indicates that the shape of the spectrum 
is accurately calculated from both Rayleigh data sets; however, slight amplitude differences do exist. The 1 kHz 
Rayleigh data still contains the shot noise contribution, resulting in a constant offset in amplitude from the CCA 
data, which is most apparent in figs. 15 and 16. The 10 kHz spectra were not calculated with the same frequency 
resolution as the 1 kHz spectra since the length of segments ( L ) must be a power of 2 and the sampling rates (f s ) 
varied by a factor of 10. The frequency “bins” for the 10 kHz spectra are nearly 2 Hz larger, resulting in slightly 
more “power” at each point in the spectrum (Fig. 16). As expected, the strongest fluctuations are found in the shear 
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layer, where the higher temperature air from the jet core is mixing with the colder co-flow air. The spectra at the 
other two locations are relatively flat. 



Figure 15. Power spectral density of temperature 
fluctuations at jet centerline. 



Frequency, Hz 


Figure 16. Power spectral density of temp- 
erature fluctuations in jet shear layer. 



Frequency, Hz 

Figure 17. Power spectral density of temperature 
fluctuations in co-flow region. 


Figures 18 and 19 show sample density and velocity PSD plots, respectively. The density PSD is shown for a 
location in the shear layer demonstrating strong density fluctuations resulting from strong temperature fluctuations. 
The velocity spectra at all locations were very similar and the centerline location was chosen as a typical case. The 
frequency scale has only been shown to 2.5 kHz to show finer detail in the low frequency range. Velocity 
fluctuations were very small due to the low speed of the jet; however, dominant fluctuations are seen at 48 and 430 
Hz, indicated by large spikes in the 10 kHz velocity PSD. The jet velocity does not actually have strong fluctuations 
at these frequencies. Rather, these spikes in the spectra can be explained by external noise sources that affect the 
frequency of the incident light. Since the velocity measurement is based on a shift in frequency from the incident 
laser frequency, fluctuations in the laser frequency become manifest in the velocity measurements. The 430 Hz 
spike was seen in previous Rayleigh measurement work in a free jet 17 using the same laser model, suggesting that 
the laser has a resonance mode around 430 Hz that is excited by the broadband jet noise. The 430 Hz modulation 
was reduced in the previous work by mounting the laser head inside an anechoic box. Modulation induced by a 
small motor used to remove speckle from the light transmitted by the fiber, often referred to as “scrambling”, may 
explain the 48 Hz spike. 
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Figure 18. Power spectral density of density 
fluctuations in jet shear layer. 



Frequency, Hz 


Figure 19. Power spectral density of velocity 
fluctuations at jet centerline. 


C. Jet profiles 

The mean temperature and velocity were calculated from the time history data, and the turbulence temperature 
and velocity fluctuations were derived from the power spectral density calculations by taking the sum of all points in 
the PSD. Figures 20 and 21 show profiles resulting from radial scans across the jet flow of the mean temperature and 
turbulent temperature fluctuations calculated from the 1 kHz CCA measurements and the 1 and 10 kHz Rayleigh 
measurements. Mean temperature values from the three data sets are found to be within +/- 3 K of each other. The 
CCA measurements provide very accurate turbulence levels for fluctuations up to 0.5 kHz; however, as the 10 kHz 
power spectra show (figs. 15-17), significant temperature fluctuations still exist beyond 0.5 kHz. Since the 1 kHz 
CCA and Rayleigh measurements do not account for higher frequency fluctuations, the resulting turbulence levels 
are lower than those calculated from 10 kHz Rayleigh data. The shot noise contribution was not removed from the 1 
kHz Rayleigh data, resulting in slightly higher fluctuation measurements than the CCA data indicates in regions 
where temperature fluctuations are small and shot noise dominates (i.e., in the jet core and in the co-flow). 




Radial position, mm 


Figure 20. Mean temperature profile at an axial 
station 32 mm downstream of the tube exit. 


Figure 21. Temperature fluctuation profile at an 
axial station 32 mm downstream of the tube exit. 


Figure 22 shows profiles of the mean velocity measured by the pitot probe and calculated from the 1 and 10 kHz 
Rayleigh data. As mentioned in section IV. A., the relative error in the velocity measurements is quite high since the 
flow velocity is so low. The error in velocity is about +/- 5 m/s while the maximum velocity is only 9 m/s. The 
results show that it is feasible to obtain velocity measurements with the presented Rayleigh technique; however the 
measurements are more applicable in higher velocity flows where the relative error is much lower. It appears from 
Fig. 22 that the velocity measurements on the lower side of the jet flow are much noisier than those obtained above 
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the jet centerline. This may be due to differences in the amount of stray light entering the optical fiber as the jet is 
translated. The amount of stray light was found to be greater for negative probe locations than for positive ones. The 
stray light level has been evaluated and accounted for in the density calculations from PMT 5. The amount of light is 
not accounted for in the least squares analysis of PMT 1-4 data since it did not seem to have a significant affect on 
the results during the calibration study. In the future, more precautions for eliminating stray light will be taken to 
improve the measurements. Figure 23 shows the turbulent velocity fluctuations calculated from the 1 and 10 kHz 
Rayleigh measurements. The actual velocity fluctuations in this jet are unknown since hotwire velocity 
measurements were not available for comparison. The turbulence levels shown here are higher than expected since 
the maximum fluctuations in jets are typically 15-20% of the difference between the centerline velocity and ambient 
or co-flow velocity. For this flow, maximum fluctuations should be on the order of 1 m/s. Sources of error in the 
velocity turbulence levels calculated from the Rayleigh data may include fluctuations induced by laser frequency 
modulation, shot noise contribution, and noise from stray laser light. 




Figure 22. Mean velocity profile at an axial station Figure 23. Velocity fluctuation profile at an axial 

32 mm downstream of the tube exit. station 32 mm downstream of the tube exit. 

V. Conclusion 

A technique for obtaining dynamic gas temperature, velocity, and density measurements using molecular 
Rayleigh scattering was described. Density was determined from an overall intensity measurement of the scattered 
light, while temperature and velocity were determined by analyzing the scattered light with a Fabry-Perot 
interferometer. The signals from five photomultiplier tubes were simultaneously recorded using photon counting 
electronics operating at 1 and 10 kHz sampling rates with 30 second recording periods. An uncertainty analysis was 
presented that demonstrated uncertainties in instantaneous temperature, velocity, and density measurements of about 
17 K, 27 m/s, and 4%, respectively, for a 10 kHz sampling rate. Derived statistical quantities, such as mean square 
fluctuations and power spectra, can be obtained with high accuracy if inherent noise contributions are eliminated, 
either experimentally or numerically in the data processing. Measurements using the Rayleigh technique were 
demonstrated in a low speed heated air flow and were verified by comparing results with CCA and pitot probe 
measurements. 

In future work, a higher power laser will be used to increase photon counts, thereby improving measurement 
uncertainty. Measurements will be made in a high speed heated flow where velocity fluctuations are significant, and 
fluctuations will be verified by hotwire velocity measurements. The presented Rayleigh scattering technique will 
eventually be used in aeroacoustics research, where sound pressure fluctuation (microphone) measurements will be 
acquired simultaneously with Rayleigh measurements allowing correlation between flow property fluctuations and 
noise generation. Determining sources of jet noise will help engineers to design quieter, more efficient aircrafts. 
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